Universal fluctuations in the relaxation of structural glasses 
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The presence of strong local fluctuations - dynamical heterogeneities - has been observed near the 
glass transitions of a wide variety of materials. Here we explore the possible presence of universality 
in those fluctuations. We compare the statistical properties of fluctuations obtained from numerical 
simulations of four different glass-forming systems: two polymer systems and two particle systems. 
We flnd strong evidence for universality, both in the qualitative behavior of the fluctuations and in 
the remarkable agreement of the scaling functions describing them. 
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Dynamical heterogeneities, i.e. strong fluctuations as- 
sociated with nanometer-scale regions of molecules re- 
arranging at very difTerent rates compared with the 
bulk [1, 2], have been observed in a wide variety of glass- 
forming systems, from small molecules, to polymers, to 
network glasses, to colloidal glasses [1, 2, 3, 4]. A de- 
tailed theoretical explanation for those strong fluctua- 
tions is not yet available, although several ideas have 
been proposed [5, 6, 7, 8, 9]. A growing experimental 
and numerical literature has been uncovering various as- 
pects of dynamical hctcrogentities [1. 2, 3, 4. 10, 11, 12, 
13, 14, 15, 16; 17, 18], but the question of universality 
-to what degree these fluctuations behave in the same 
way in different kinds of glassy systems- remains open. 

Another common feature of glass-forming materials is 
physical aging [19]: for temperatures below the glass 
transition, the material falls out of equilibrium, and 
quantities probing the system at two times, a "waiting 
time" t^ and a "final time" t, with t > tm, depend on both 
times tu, and t, and not just on their difference t — t^^. 
In particular, probability distributions and spatial cor- 
relations describing dynamical heterogeneities also show 
aging [10]. Their time dependences [11, 12, 13, 14] 
display scaling as a function of the two-time correla- 



tion C, 



global 



, t^ ) 



E^iexp(iq.(rj(t) - rj(t„))), 



as predicted by a theoretical framework based on the 
presence of local fluctuations in the age of the sam- 
ple [9, 20, 21, 22, 23]. 

In this work we examine the question of universality 
in the fluctuation behavior of structural glasses. We per- 
form numerical simulations in four different models of 
glass-forming systems in the aging regime, and take ad- 
vantage of the presence of scaling to quantitatively com- 
pare the properties of fluctuations in all of them. As a 
result of those comparisons, we find strong evidence that 
dynamical heterogeneities do indeed exhibit universal be- 
havior. 

We consider a system of polymers (labeled Ijp) with 
Lennard Jones (LJ) interactions between the monomers, 
a system of polymers (labeled Wp) with purely repulsive 



Weeks-Chandler- Andersen (WCA) interactions [24, 25], 
a system of particles (labeled Ijm) with LJ interactions, 
and a system of particles (labeled w„i) with WCA in- 
teractions. In our polymer models, nearest-neighbor 
monomers along a chain are connected by a FENE an- 
harmonic spring potential, and we use the Nose-Hoover 
method to simulate at a constant pressure and temper- 
ature [24]. Each system, composed of 800 chains of 10 
monomers each, is equilibrated at a high temperature of 
Ti = 5.0 and then it is instantaneously quenched to a final 
temperature Tf . All temperatures are measured in units 
of the energy scale e of the LJ or WCA potential. The 
time of the quench is taken as the origin of times t = 0. 
After the quench the systems are allowed to evolve for 
10^ LJ time units. Tf is chosen low enough that the sys- 
tems keep aging during the whole low temperature part 
of the simulation: Tf = 0.4 for Wp, and Tf = 0.6 for Ijp. 
Our results are an average over 100 (resp. 800) indepen- 
dent runs for the Ijp (resp. Wp) system. The simulations 
in the particle systems are as described in Refs. [12, 14]. 

We present results for the probability distributions of 
observables which probe local fluctuations in small re- 
gions of the system: the local coarse grained two-time 
correlation function [12, 14] Cr(t, tto) and the particle 
displacements along one direction Axj{t,tw) = Xj{t) — 
Xj{tiu) [4, 12, 15]. In order to probe the spatial correla- 
tions of the fluctuations, we also consider the generalized 
dynamic susceptibility Xi = J d'^^ 54(r,t,t«;) [13, 16, 17, 
18], where (?4(r, t, t^u) is a 4-point (2-time, 2-position) cor- 
relation function [13, 16]. 

The top panel of Fig. 1 shows the probability dis- 
tributions p(Cr) of the local two-time correlation for 
Cgiobai(i,t«;) = 0.2,0.4,0.6, for WCA polymers. As in 
particle systems [12, 14], the probability distributions 
are approximately independent of tw, for a fixed value 
of Cgiobai(i, ^tii)- I'^ the bottom panel of Fig. 1 we plot 
the rescaled probability distributions crcp(Cr) versus the 
normalized fiuctuation (Cr — CgiobaO/cc in the one-point 
two-time correlator for Cgiobai(i, t^) = 0.5, and find that 
the results are approximately the same for all four sys- 
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FIG. 1: Probability distribution p{Cr), for coarse grain- 
ing regions containing on average 6.6 particles. Top panel: 
p(Cr) in the aging regime of a WCA polymer glass when 
Cglobal(ii = 0.2,0.4,0.6. Bottom panel: Rescaled prob- 
ability distributions ocp{Ct) as functions of the normal- 
ized fluctuation 
for 
m 

except t 



(Cr - C„iobal)/o-C for Cglobal = 0.5, 



the four systems, with q = 7.24,7.37,7.20,7.20 for 
p, ljp,Wm, Ijm respectively. Here t^^ — 30.2, 3020 in all cases 



= 30.2 only for WCA particles. 



terns: Ijm, Wm, Ijp ^nd Wp. 

The top panel of Fig. 2 shows the probability distri- 
butions p{Ax) of particle displacements for the four sys- 
tems, for fixed Cgiobai(i, iuj) = 0.5 and = 30.2,3020. 
The distributions look similar to those determined by 
confocal microscopy in a colloid with repulsive interac- 
tions [4] but rather different from those found in an at- 
tractive colloid [15]. All the distributions collapse to- 
gether to a common Gaussian shape for smaller Ax, but 
depart from the Gaussian for larger Ax. For the shorter 
waiting times, ~ 30.2, the tails of the distributions in 
the four systems arc very close to each other, but they sig- 
nificantly separate for longer tw In all systems the tails 
of the distributions become wider as increases; this 
is shown in more detail for the Wp system in the middle 
panel of Fig. 2. In all four systems, the tails can be fit 
in the region |Aa;| > 0.5 by a nonlinear exponential form 
p(Ax) « A^exp(— |Ax/a|'^). However, the bottom panel 
of Fig. 2 shows that for WCA polymers, the exponent 
increases at long t^, while for LJ particles, it decreases 
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FIG. 2: Top panel: Non-universality of the time evolution of 
the tails of p{Ax). Plot of p{Ax) for Cgiobal(^i = 0.5 for 
the four systems, with t^n ~ 30.2, 3020 (except t^ — 30.2 only 
for WCA particles). Plotted with a logarithmic vertical axis 
to emphasize the tails of the distributions. Gaussian fits to all 
seven curves are also shown (black full lines); all fits collapse 
with each other. Middle panel: p{Ax) with Cgiobal(^i ^m) = 
0.5 and t^ = 30.2, 302, 3020, 30200; for a WCA polymer glass. 
Bottom panel: Evolution with waiting time of the (3 exponent 
describing the tails of p{Ax), for a polymer system and a 
small molecule system, at Cgiobai(*i ^m) = 0.1,0.3,0.5. 



at long tjjj. According to [15] the tails of the distribution 
are representative of particles which have escaped their 
cages, whereas the peak of the distributions are due to 
the particles which vibrate in place. This suggests that 
the difference in the evolution of the tails could be due to 
the different diffusive behavior between small molecules 



3 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.2 0.4 0.6 0.8 1 

1-C(t,tJ 1-C 



FIG. 3: Top panel: Xii't^'tw) as a function of t — t^i, with 
constant = 100, ■ • • , 30200, for WCA polymers (data from 
5000 independent simulation runs). Middle panel: Rescaled 
X4, plotted as a function of 1 — C, with constant — 
30.2, • ■ ■ , 30200, for WCA polymers. A very good collapse 
is observed. Bottom panel: Rescaled X4i plotted as a func- 
tion of 1 — C, with various waiting times, for three of the four 
systems. 



and polymers. 

We now consider the spatial correlations of dynam- 
ical fluctuations. The top panel of Fig. 3 shows that 
Xiiti tw) has a peak as a function of t — t^, and that this 
peak's height and position grow with . Similar behav- 
iors are observed in experiments in granular systems [18] 
as the area fraction is increased, and in numerical sim- 
ulations of supercooled liquids as the temperature is re- 
duced [16, 17]. In the middle panel, wc plot the ratio 



FIG. 4: Top panel: Rescaled S4{q,t,tw), as a function of 
the scaling variable x = ^4,{t,tui)\q\- Data are shown for 165 
time pairs (t, tm) for LJ particles and 338 time pairs for WCA 
polymers. Inset: ratios of the average squared error divided 
by its minimum value e^jj^ for Ijm and Wp systems. The curve 
labeled Wp + Ijm shows the average of the two ratios. Middle 
panel: £,4{t,tw) as a function of t — t^, with constant t^ ~ 
100, ■ ■ ■ , 30200, for WCA polymers. The value of ^4 at the 
plateau increases with t^. Bottom panel: Rescaled ^4,{t,tw), 
as a function of 1 — C The upper plot is for LJ particles with 
tjn = 13.2, ■ ■ ■ , 832 and the lower plot corresponds to WCA 
polymers with constant t^ — 100, • ■ ■ , 30200. 



X4/X4(C=l/e) as a function of 1 - C = 1 - Cgiobai(*, i™), 
for WCA polymers; all the curves collapse into a single 
master curve. Both behaviors are identical to those found 
in LJ and WCA particles [13, 14]. In the bottom panel 
of Fig. 3, the same ratio is plotted for L J particles, WCA 
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particles and WCA polymers. The scaling function is 
similar but apparantly not identical in the polymer and 
small molecule cases. 

By Fourier transforming the correlation function 
54 (r, t^tw), we obtain the 4-point dynamic structure fac- 
tor S'4(q, t,t„) [13, 16, 17]. We fit its small q behav- 
ior (g < 1.9) with an empirical scaling form: 5*4 (q) = 
{Si - K) f{UH) + K, with f{x) = 1/(1 + x'^Y, for the 
Wp and Ijm systems. The dynamic correlation length 
i,i{t,tyj) and the parameters Sl{t,t^) and K(t,tw) are 
extracted from each fit. The unknown scaling function 
f{x) is the same for all the fits. It has sometimes been 
assumed [16] that it has an Ornstein-Zernicke (OZ) form, 
corresponding to /i = 1. More recently, however, it has 
been shown [26] that the function Xq, which is believed 
to have a similar behavior to 5*4 (q), does not have an OZ 
form, and indeed has an asymptotic behavior Xq ~ 
for large q. Indeed, forms other than OZ are better at 
describing the data in the supercooled regime of various 
glassy models [17], and they provide significantly better 
fits to the aging regime data presented here. By mini- 
mizing the average square error over the whole set of fits, 
we determine ^ ~ 1.8 [27]. Although our fits are per- 
formed only for moderate to small values of x = ^4(7, the 
asymptotic behavior for large argument, f{x) ~ a;"'^'^, is 
not far from the prediction of [26] . 

The top panel of Fig. 4 shows that all the data sets 
agree rather well with the scaling function. The middle 
panel of Fig. 4 shows the extracted correlation lengths 
vs. t — tyj for different waiting times, for WCA polymers. 
We see an initial increase in the correlation length in all 
the curves, but when the time difference gets larger, the 
correlation length either remains constant or it decreases 
slightly [28]. As in LJ particles [13], the plateau value 
of ^4 is a growing function of t.^, . The bottom panel of 
Fig. 4 shows plots of the ratio £,i/£,i{C = 1/e) against 
1 — C, for the Wp and Ijm systems. In both cases we find 
that there is a moderately good collapse, and it appears 
that ^4 goes to a nonzero constant when (1 — C) 1, 
(i.e. {t — tw)/tu! — > 00). However, the scaling function for 
monomers grows with 1 — C for 1 — C > 0.5, while for the 
polymers it appears to become approximately constant. 

In summary, we have explored the possible presence 
of universality in the fluctuations of relaxing structural 
glasses, by comparing simulation data from two poly- 
mer models and two particle models. As we pointed 
out before, some differences are observed: (a) the scal- 
ing plots for Xi/Xii^ = ^-iid ii/^iiC ~ 1/e) versus 
1 — C show some small discrepancies between polymers 
and particles; and (b) the behavior of the tails of the dis- 
tributions of displacements, p{Ax(t,tw)), show a quali- 
tatively different evolution with ty^ in the small molecule 
and polymer cases (which might be due to their different 
diffusive behaviors). However, there is remarkable sim- 
ilarity in the behavior of fluctuations in all the systems 
considered, and we find the evidence for universality to 



be very strong. In particular, the following aspects of 
the fluctuations appear to be universal: (i) the fact that 
the probability distributions p{Cr{t^tw)) approximately 
collapse for different waiting times when C(t,tw) is 
held constant; (ii) the very peculiar shape of the scal- 
ing function p((Cr - Cgiobai )/o-c) = crcp{Cr); (iii) the 
qualitative behavior of the 4-point density susceptibility 
X4(<,tu,) and the dynamic correlation length S,4:{t,tw) as 
functions of and t ~ t^; and (iv) the fact that the 
rescaled quantities X4/x4(C' = 1/e) and ^4/^4 (C = 1/e) 
plotted versus 1 — C approximately collapse for different 
waiting times 
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